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Abstract 

We extend the powerful frequency-domain analysis of femtosecond two-photon absorption to the 
intermediate-field regime, which involves both two- and four-photon transitions. Consequently, 
we find a broad family of shaped pulses that enhance the absorption over the transform-limited 
pulse. It includes any spectral phase that is anti-symmetric around half the transition frequency. 
The spectrum is asymmetric around it. The theoretical framework and results for Na are verified 
experimentally. This work opens the door for rational femtosecond coherent control in a regime of 
considerable absorption yields. 
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Following from their coherence over a broad spectrum, shaped femtosecond pulses allow 
to coherently control quantum dynamics in ways that other means cannot [H [21 El HI E]- 
Among the processes, over which such control has been most effective, are the important 
multiphoton absorption processes [HI El El El El HOj HH [121 H31 H31 UHl [161 HZl HE! H9] - 
The control is achieved by shaping the pulse to manipulate the interferences between vari- 
ous initial-to-final state-to-state multiphoton pathways. Constructive interferences enhance 
the absorption, while destructive interferences attenuate it. Thus, to fully utilize coher- 
ent control, the pulse shaping should ideally be based on identifying first the pathways 
and their interference mechanism. When such identification is not feasible, a practical par- 
tial solution is to employ automatic "black-box" experimental optimization of the pulse 
shape using learning algorithms [20]. For multiphoton processes, the past control studies 
[HI [61 [71 El IS1 HOI HH US! H3l H^l US] have shown that this ideal line of action is feasible 
and very powerful once the photo-excitation picture is available in the frequency domain. 
However, so far the frequency domain has been exploited only in the weak-field regime 
[3 El El El O HOI [HI H21 H31 [HI HH] associated with low absorption yields (typically up to 
~0.1% population transfer). There, the N-photon absorption is described by time- dependent 
perturbation theory of the lowest non-vanishing order (the N th order) and, thus, can be 
transformed to the frequency domain. 

The present work extends the powerful frequency-domain picture to a regime of sizeable 
absorption yields, exceeding the weak-field yields by two orders of magnitude (and more). 
In this regime the interfering pathways of N-photon absorption are the weak-field pathways 
of N absorbed photons as well as additional pathways of M absorbed photons and M— N 
emitted photons (M>N). The corresponding picture is obtained by extending the perturba- 
tive analysis to include a finite number of non-vanishing orders beyond the lowest one. We 
refer to this regime as the intermediate-field regime. It is distinguished from the strong-field 
regime where no perturbative description is valid, which is the one all the past multiphoton 
control studies beyond the weak-field regime have focused on [161 El EIE1 HE] • Specifically, 
based on 4 i/l -order perturbation theory, we develop here an intermediate-field frequency- 
domain theoretical framework to the process of atomic femtosecond two-photon absorption, 
which is non-resonant in the weak- field regime (of 2 nd order). As a 4 i?t -order approach, it 
involves both the 2 nd and 4 th perturbative orders associated, respectively, with two- and 
four-photon pathways. Consequently, we find an extensive family of shaped pulses that in 
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this regime enhance the absorption as compared to the unshaped transform-limited (TL) 
pulse. This is impossible in the weak-field regime [Bj. The theoretical framework and results 
are verified experimentally. The approach developed here is general and can be employed in 
other cases. 

The physical model system of the study is the sodium (Na) atom (see Fig. Q. Theo- 
retically, we consider an atomic two-photon absorption process from an initial ground state 
\g) to a final excited state |/), which are coupled via a manifold of intermediate states \n) 
having the proper symmetry. The spectrum of the pulse is such that all the |<7)-|ji) and 
\f)-\n) couplings are off-resonant, i.e., the spectral amplitude at all the corresponding tran- 
sition frequencies is zero: |.E(u; sn )| = \E(ujf n )\ = 0, except for the \f)-\n r ) resonant coupling 
for which \E(uf nr )\ ^ 0. In the present intermediate-field regime, the final (complex) am- 
plitude Af of state |/), following irradiation with a (shaped) temporal electric field e(t), 
can be validly described by 4 t/l -order time-dependent perturbation theory. So, in general, it 
includes non-vanishing contributions from both the 2 nd and 4 th perturbative orders: 

A, = Af + Af. (1) 

Within the frequency-domain framework, the spectral field E(u) = |J5(o;)| exp [i<&(a;)] is 
given as the Fourier transform of e(t), with |-E(cu)| and $(u;) being the spectral amplitude 
and phase at frequency u. For the TL pulse: $(a;) = 0. We also introduce the normalized 
spectral field E(cj) = E(uS)/ \E \ representing the pulse shape, with \E \ being the maximal 
spectral amplitude. As shown for the weak-field regime [6J , the 2 nd -order term A^p interferes 
all the \g)-\f) pathways of two absorbed photons (see examples in Fig. [T]) and is given by 

Af—^I^W (2) 

/oo ^ 
E(uj)E(Vl - u)duj , (3) 
-oo 

where Uf g and n"j g are the \g)-\f) transition frequency and effective non-resonant two-photon 

(2) 

coupling. The TL pulse induces fully constructive interferences within A^ and, thus, the 
and maximal weak- field non-resonant two-photon absorption [6j. The 4 th - 



maximal 



A f 



order term A^ is more complicated than and we have analytically calculated it to be 
given by 



4 (4) - 1 IF I 4 



/OO -I 
d5-A^\to fg -5)A^ R \5) 



(4) 
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where A^ 2 \il) is defined above, and 

A {R) {AQ) = A (non - rcsR \AQ) + A {rcsR \AQ) (5) 

/oo 
E(co + AQ)E*(co)du (6) 
-oo 



/+oo 1 ^ ^ 

dd'-E(u fnr + AQ- 5')E*(uj fnr - 5') 
-oo " 



(7) 



These equations reflect the fact that A^ interferes all the four- photon pathways from \g) 
to \ f) of any combination of three absorbed photons and one emitted photon. 

Figure [T] shows several representative four-photon pathways. Each pathway can be divided 
into two parts: a non-resonant absorption of two photons with a frequency sum of Q = ujf g —5 
and a Raman transition of two photons with a frequency difference of AQ = 5. Their border 
line is detuned by 5 from either |/) or \g) according to whether, respectively, the two- 
photon absorption part or the Raman part comes first (see Fig. [TJ. The on-resonant (5 = 0) 
and near-resonant (5 ^ 0) pathways are interfered separately in Eq. (|4]). The Cauchy's 
principle value operator p excludes the on-resonant pathways from the second term. The 
integration over the pathways is expressed using two parameterized amplitudes, A^iVL) 
and A^ R \An), that originate from the different parts of the four-photon pathways. A^ 2 \il) 
interferes all the two-photon absorption pathways of frequency £1. A^ R \AQ) interferes all the 
Raman pathways of frequency AQ and includes two components. The first is A^ non ~ resRS) (AO) 
interfering all the Raman pathways that are non- resonant, with n 2 gg and fij^ being the \g)- 
\g) and |/)-|/) effective non-resonant Raman couplings due to the non-resonantly coupled 
states \n). The second is A^ esR \AVL) interfering all the Raman pathways that are resonance- 
mediated via \n r ), with /ij nr being the \ f)-\n r ) dipole matrix element. 

For a given pulse shape E(u), a non-zero A^ is proportional to the maximal spectral 
intensity Jo = \E \ 2 while a non-zero A^ is proportional to Iq = \E \ 4 . Different intensity 
Io corresponds to a different temporal peak intensity of the TL pulse (/tl)- In this work, 
for a set of intensities Jo, the |/)'s final population Pf = \Af\ 2 , which reflects the degree of 
two-photon absorption, is controlled via the pulse shape E(u). The spectral phases $(u;) 



are the corresponding control knobs, with the spectrum 



E(io) 



unchanged. The latter is 



chosen such that the central frequency u is detuned from ojf g /2 and E(ujf g /2) 



discussed below, this corresponds to the interesting case of having A^ negligible relative 
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to Ajr in the weak-field limit and comparable to in the upper intermediate- field limit. 
Their interplay enables the absorption enhancement beyond the TL absorption. 

The Na system [22] includes the 3s ground state as \g), the 4s state as |/), the manifold of 
p states as the \n) manifold, and the 7p state as \n r ). The transition frequencies Uf g = tt>4 Sj 3 S 
and ujf rir = uj-jp^s correspond, respectively, to two 777-nm photons and one 781. 2- nm photon. 
The sodium is irradiated with phase-shaped linearly-polarized femtosecond pulses having an 
intensity spectrum centered around 779.5 nm with 5 nm bandwidth (~180 fs TL duration). 
Experimentally, a sodium vapor in a heated cell is irradiated with such pulses of variable 
energy after they undergo shaping in an optical set-up with a pixelated liquid-crystal spatial 
light phase modulator [21] . Upon focusing, the TL temporal peak intensity at the peak 
of the spatial beam profile /te™ fiIe ~P eak ) ran ges from 10 9 to 5xl0 10 W/cm 2 . Following the 
interaction with a pulse, the population excited to the 4s state undergoes cascaded decay 
to the 3s state via the 3p state. The corresponding 3p-3s fluorescence serves as the relative 
measure for the final 4s population Pf = P 4s . It is measured using a spectrometer coupled 
to a time-gated camera system. The measured signal results from an integration over the 
spatial beam profile. 

We validate our intermediate-field frequency- domain perturbative description by com- 
paring corresponding results to exact non-perturbative results obtained by the numerical 
propagation of the time-dependent Schrodinger equation (using the Runga-Kutta method). 
The considered manifold of p-states is from 3p to 8p (22]- The validity of the non-perturbative 
calculations is confirmed first by a comparison to experiment. The test case is the set of 
shaped pulses having a 7r-step spectral phase pattern. Each such pattern is characterized by 
the step position to s tep, with $(cj < cu step ) = —n/2 and $(u; > uj step ) = ir/2. Figures |2^a)-(c) 
show examples of experimental (circles) and non-perturbative theoretical (solid lines) results 
of 

P4.S db function of 00 step for different pulse energies, i.e., different j^P rofile -P eak ) _ E ac h of the 
traces is normalized by of the TL excitation. The weak-field trace [H] is shown in Fig.[2^a). 
The theoretical traces account for the experimental integration over the beam profile: Each 
of them results from a set of calculations conducted each with a different (single) value of Iq. 
As can be seen, the agreement between the experimental and non-perturbative theoretical 
results is excellent, confirming these calculations' accuracy. Figs. |2^d)-(g) compare, for dif- 
ferent (single- valued) intensities J , the theoretical non-perturbative results (squares) with 
perturbative results calculated numerically using Eqs. (|l])-([7]) (solid lines). Shown are exam- 
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pies out of the full set of results. As can be seen, the perturbative results reproduce the non- 
perturbative ones up to the Iq corresponding to a TL peak intensity of Jtl=2.7x 10 10 W/cm 2 
[Fig. pm]. This is the intensity limit of the intermediate-field regime for the present Na 
excitation, up to which no perturbative order beyond the 4 th one is needed to be included. 
The weak-field regime of A^^A^J extends here up to Jtl~10 9 W/cm 2 [Fig. 2^d)]. 

As shown in Figs. |2jd)-(f), the 7r-trace shape significantly changes as Iq increases. The 
prominent feature reflecting the deviation from the weak-field regime is the absorption en- 
hancement beyond the TL absorption when uj step is around u)4 S ^ s /2 (777 nm) or around 
the central spectral frequency Uq (780 nm). The enhancement increases as Iq increases. 
Here, using the intermediate-field description, we analyze the enhancement mechanism 
when uj s tep = ^4s,3s/2. In the weak-field regime [Fig. [2[d)] the corresponding absorption 
is equal to the (maximal) TL absorption, since such a phase pattern also leads to fully con- 
structive interferences within A^J [6]. In the present intermediate-field limit [Fig. 2^f)] 
it is about twice the TL absorption: With the Iq of Jtl=2.7x10 10 W/cm 2 the calcu- 
lated values are {AfJ = 0.35, = -0.15 - 0.10* : P 4s = 0.05} for the TL pulse and 
{A^ ] s = 0.35, A$ = -0.09 - O.Ui : P As = 0.09} for the pulse of u step = uj AS)3s /2. So, as seen 
from these values, the lower TL absorption results mainly from a stronger attenuation of 
the equal real A A ] by 9ft [A^]. 

The difference in 3ft[^4^] for these two pulses originates from the jp-integral in Eq. (4), 
which interferes all the near-resonant four-photon pathways of 5 7^ with the domination 
of small \5\ (due to the 1/5 weighting). Actually, it predominantly originates from the 
integrand factor A^^uo^^g — 5) interfering all the two-photon absorption parts (of these 
four-photon pathways) of transition frequency uj As ,zs — 5. For the TL pulse, due to fully 
constructive interferences, A^ (to As ,3s — 5) is maximized (real and positive value) for any 5. 
So, with a detuned spectrum of u Q ^ u iSi3s /2, it changes monotonically with 8 around 5 = 0. 
A^ (lu^s^s — \6\) > v4( 2 )(cg>4 Sj3s + \5\) with a red spectral detuning, and vice versa with a 
blue detuning. On the other hand, for the shaped pulse of uj s te P =uJ ASl 3 S /2, A^ 2 \uj iSt ss — 5) 
is maximized only for 5 = and gradually reduces for small non-zero |<5| with comparable 
magnitude for ±\5\ (true for both red and blue spectral detuning). Consequently, considering 



(4)1 
4s J 



for the TL 



also the 1/5 integrand factor, the jp-integral results in a higher value of 9ftL4 
pulse. For both pulses, the sign of 3ft [A^] relative to A® is determined by the pulse spectrum 
and by the Raman couplings of the four-photon pathways (/i| s3s , /x| s 4s and |/U 4S)7 p| 2 ). Here, 
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the present red-detuned spectrum (uj of 779.5 nm) leads to and 9ft [A^J] that are of 
opposite signs. Shifting the spectrum to be comparably blue detuned from u is ^ s /2 (u of 
774.5 nm) changes their signs to be the same. Then, the intermediate-field absorption of the 
TL pulse is the higher one. A non-detuned spectrum (ujq of 777 nm) leads to a negligible 
] and thus to an equal absorption for both pulses, in agreement with the strong- 
field time-domain study by Dudovich et al. [16J. This spectral dependence is the subject 
of future publication. For the current red-detuning, the above explanation also allows to 
understand why the enhancement over the TL absorption increases as Jo increases: Since 
SftLA^] /A4J = KI with higher value of K for the TL pulse as compared to the shaped 
pulse of uj ste p=uj is ^ s /2, a given increase of Iq results in a higher increase in the attenuation 
of A^ by A±g for the TL pulse. 

Similar to the 7r-step phase pattern of u step =uj4 :St3s /2, all the spectral phase patterns 
that are anti-symmetric around U4 St 3 S /2, i.e., $(u;) = — $(co>4 Si 3 S — u), lead to the maxi- 
mal weak-field two-photon absorption for, and only for, a transition frequency of uJ4 S .3s [§]• 
So, based on the above analysis, in the intermediate-field regime the corresponding shaped 
pulses are expected to enhance the two-photon absorption (to the 4s state) beyond the TL 
absorption. The mechanism is the enhanced attenuation of A^J by A^j for the TL pulse. 
Figure [3] confirms that indeed this is the case. It shows theoretical [Fig. ^a)] and exper- 
imental [Fig. |3](b)] results for 5000 different pulse shapes with such anti-symmetric phase 
patterns that we have randomized. The results are presented as histograms (distributions) 
showing the fraction of pulses inducing different values of P4 S /PiJ L \ i-e., different absorp- 
tion enhancement over to the TL excitation. The theoretical results are perturbative and 
have been calculated numerically using Eqs. ([T])-([7]), with each histogram corresponding to 
a different (single- valued) intensity Iq. The (spatially-integrated) experimental histograms 
correspond to different pulse energies. The weak-field histogram is located at the value of 1. 

As can be seen, the distribution shifts to higher enhancement values and gets broader as 



I increases. Both effects originate from the J -dependence of 9ft [A^] /A^J = KI , with K 
depending only on the pulse shape. Similar to the uj s te P =0J4 LS fi S /2 case, the shift to higher 
enhancement values originates from the higher K value for the TL pulse as compared to the 
shaped pulses with the anti-symmetric phase patterns. The distribution broadening origi- 
nates from the variation in the K value from one pulse shape to the other. The experimental 
histograms are generally broader than the theoretical ones due to the integration over the 



I (2) 
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spatial beam profile and due to the experimental noise. The experimental signal-to-noise ra- 
tio is directly reflected in the width of the weak- field histogram [Fig. [3|b)-dashed line]. Our 
most-enhancing random anti-symmetric pattern is calculated at Io of 7tl=2.7x10 10 W/cm 2 
to induce P4 S =0.21, corresponding to 4.2 enhancement over the TL absorption. Last, it 
worth mentioning that excluding the resonantly-coupled 7p state from our calculations gen- 
erally leads to about 25% reduction in the enhancement factors over the TL absorption at 
J TL =2.7xl0 10 W/cm 2 . 

In summary, we have extended the powerful frequency- domain analysis of femtosecond 
two-photon absorption to the intermediate-field regime of 2 nd and 4 th perturbative orders. 
This opens the door for rational femtosecond coherent control beyond the weak-field regime, 
in a regime of considerable absorption yields. The extended frequency- domain description 
has enabled us to identify a broad family of shaped pulses that enhance the absorption 
beyond the TL absorption. The mechanism is the coherent interplay between the different 
perturbative orders. We expect this work to serve as a basis for intermediate-field extensions 
to molecules and other multiphoton processes as well as for devising effective control strate- 
gies when several multiphoton channels are involved. The results are also of applicative 
importance for nonlinear spectroscopy and microscopy. 
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FIG. 1: (Color online) Excitation scheme of the intermediate- field two-photon absorption in Na 
(not to scale), with representative two-photon and four-photon pathways from \g) = 3s to |/) = 4s. 
Each of the four-photon pathways is either on- or near-resonance with 4s or 3s. Its Raman part is 
non-resonant due to the np states (n^=7) or resonance-mediated via 7p. 
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FIG. 2: Results for the Na two-photon absorption induced by the shaped pulses with a it spec- 
tral phase step as a function of the step position. Panels (a)-(c): Experimental (circles) and 
non-perturbative theoretical (solid lines) results (spatially integrated) for different pulse energies 
corresponding to a TL intensity at the peak of the spatial beam profile j^ ofile "P eak ) Q f ( a ) iq9^ ^ 
2.3xl0 10 , and (c) 3.4xl0 10 W/cm 2 . Panels (d)-(g): Non-perturbative (squares) and perturbative 
(solid lines) theoretical results for different (single-valued) spectral intensities 1$ (see text) corre- 
sponding to a TL intensity J TL of (d) 10 9 , (e) 1.6xl0 10 , (f) 2.7xl0 10 , and (g) 4xl0 10 W/cm 2 . The 
perturbative calculations include 2 nd and 4 th orders. 
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FIG. 3: Histograms of the Na two-photon absorption for a set of 5000 different pulses with random 
spectral phase patterns that are all anti-symmetric around U4 S ^ s /2. They show the fraction of 
pulses inducing different absorption enhancement over the TL pulse, (a) Theoretical perturbative 
results (of 2 nd and 4 th orders) for different (single-valued) spectral intensities Iq corresponding 
to different TL intensities Itl- (b) Experimental results (spatially integrated) for different pulse 
energies corresponding to different TL intensities at the peak of the beam profile j&-° me -P ca -k) ^ 
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